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ABSTRACT 

We present new long-slit William Herschel Telescope spectroscopic 
observations of the lens galaxy Gl associated with the double- imaged 
QSO 0957+56lA,B. The obtained central stellar velocity dispersion, 
oi = 310 ±20 km s" 1 , is in reasonable agreement with other measurements of this 
dynamical parameter. Using all updated measurements of the stellar velocity 
dispersion in the internal region of the galaxy (at angular separations < 1'.'5) and 
a simple isotropic model, we discuss the mass of a possible central massive dark 
object. It is found that the data of Falco et al. (1997) suggest the existence of 
an extremely massive object of (0.5-2.1) x 10 10 hr l M & (80% confidence level), 
whereas the inclusion of very recent data (Tonry & Franx 1998, and this paper) 
substantially changes the results: the compact central mass must be < 6xl0 9 
h~ x M Q at the 90% confidence level. We note that, taking into account all the 
available dynamical data, a compact nucleus with a mass of 10 9 hr x M Q (best 
fit) cannot be ruled out. 

Subject headings: dark matter - galaxies: elliptical and lenticular, cD - galaxies: 
structure - galaxies: kinematics and dynamics - gravitational lensing - quasars: 
individual (0957+561) 
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1. Introduction 

Gravitational lenses have turned out to be ideal sites for obtaining estimates of several 
cosmological quantities, such as the content of the universe (see, for example, Im et al. 

1997) or the Hubble constant. Refsdal (1964a,b) first demonstrated that a determination of 
H could be obtained from the time delay between two components of a multiply-imaged 
quasar and an estimate of the distribution of mass in the deflector. This method is based 
on purely geometrical grounds (within the framework of the general theory of relativity), 
and presents decisive advantages with respect to conventional astronomical techniques (see, 
for example, Kundic et al. 1997). In particular, it avoids uncertainties in the construction 
of the cosmological distance ladder, and allows the expansion rate from high-redshift — i.e., 
cosmologically distant — objects to be obtained. 

Q0957+561A,B was the first gravitational lens system to be discovered (Walsh, 
Carswell, & Weymann 1979), and has been monitored since then (see, for example, Schild & 
Thompson 1995; Oscoz et al. 1996). However, despite the several observational campaigns 
carried out, the value for the time delay was not well determined until very recently 
(Kundic et al. 1997; Oscoz et al. 1997). Nowadays, a time delay close to 420 days is widely 
accepted. Furthermore, it is known that the main gravitational effect in Q0957+561A,B is 
produced by two intervening objects: a giant elliptical galaxy (Brightest Cluster Galaxy, 
Gl), discovered by Stockton (1980), and a cluster of galaxies surrounding Gl. Several 
mass models, in which the form of the lens (galaxy + cluster) profile were discussed, have 
been proposed since 1979 (Young et al. 1980; Borgeest & Refsdal 1984; Greenfield et 
al. 1985; Kochanek 1991; Falco et al. 1991; Grogin & Narayan 1996a,b; Barkana et al. 

1998) . Unfortunately, the models are incomplete (the fits are unable to discriminate the 
contributions from mass associated with the galaxy and mass associated with the cluster), 
and are therefore unable to compute the Hubble constant without additional information. 
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To break this degeneracy, a measurement of either the velocity dispersion of Gl (the ID 
velocity dispersion associated with the dark objects in the halo of Gl) or the surface mass 
density of the cluster (Fischer et al. 1997) is needed. 

The first attempt to measure 07 (the ID velocity dispersion associated with the 
luminous stars in Gl) was made by Rhee (1991) who, analyzing the Mg lb triplet absorption 
lines from the spectra obtained with the 4.2-m William Herschel Telescope (WHT), inferred 
a value for the stellar velocity dispersion of Gl of 07 = 303 ± 50 km s -1 . Some years later, 
Falco et al. (1997, hereafter FSMD) have also measured the Mg lb triplet with the Keck 
Telescope, reporting an apparent decrease of 0/ with angular distance from the center of 
Gl (from 316 km s" 1 at ip w Of 1 to 262 km s" 1 at if p=* Of 5) and obtaining an averaged 
value of 279 ± 12kms _1 . Any giant elliptical could have a central black hole of 10 7 — 10 9 M Q 
(see Kormendy et al. 1998), and so the gradient has been preliminarly interpreted by the 
authors as due to a central massive dark object of mass 4x10 9 hr x M Q (H = 100 h km 
s _1 Mpc -1 ), i.e., a mass similar to the largest measured central mass in elliptical galaxies. 
By means of a more careful analysis (see §4.2), it can be shown that the decrease of 07 with 
distance to the galaxy center reported by FSMD should imply a central object with a very 
large mass of ps 10 10 hr l M Q (assuming, as do FSMD, that the gradient is only due to a 
compact nucleus). Very recently, Tonry & Franx (1998, hereafter TF) have measured (also 
with the Keck Telescope) the central stellar velocity dispersion inside an aperture of l"x 
2f3 as well as the 07 from different rows along the spectrograph slit. The slit was rotated by 
90 deg with respect to the previous orientation by FSMD, and with this new configuration, 
no gradient is seen in the stellar velocity dispersion trend. Therefore, TF do not agree 
with the conclusion of FSMD concerning the existence of a central dark object with very 
large mass. In fact, their dynamical data are in apparent agreement with the absence of a 
compact nucleus. 
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We present here a new measurement of <T/ (the spectrograph slit was rotated by about 
10 deg and 20 deg with respect to the exposures 1—3 and 4 by FSMD, respectively) of 
similar accuracy to those of FSMD and TF. In §2 we present the observations and data 
reduction. Section 3 is devoted to the derivation of 07 {ip < 1'.'5) from the spectra, and in 
§4 we use all the updated information (the data of FSMD and TF, and ours) to discuss the 
mass of a possible central massive dark object. 

2. Observations and Data Reduction 

The spectra of the lens galaxy in Q0957+561 were collected at the WHT (La Palma, 
Canary Islands, Spain) on 1997 February 4. We used the red arm of the Intermediate- 
dispersion Spectrograph and Imaging System (ISIS) with a Tek CCD camera. The 158 
line mm -1 grating, R158R, gave a dispersion of 1.21 A pixel -1 , covering a spectral range 
from 5000 to 8000 A. The seeing was 1" FWHM and the spectral resolution, derived from 
comparison lamps, was ~ 2.2 A. We centered a l"x 240" slit on Gl (but an aperture of only 
l"x 3" is used), and oriented it along P. A. ~ 190° to include image B also in the aperture 
(see Fig. 1). We took 11 exposures of 1800 s each. We also obtained three exposures of 10 
s of the calibration star AGK 2+14783. 

Data reduction was carried out using the standard software package iraf[|. The 
spectrum images were bias-subtracted and trimmed, flat-field corrected and cleaned 
of cosmic rays. Wavelength calibration and distortion correction were performed with 
polynomial fits to the Cu-Ne arc lamp spectra. Finally, we co-added all the 1800-s spectra. 

1 iraf is distributed by the National Optical Astronomy Observatories, which is operated 
by the Association of Universities for Research in Astronomy, Inc. (AURA) under 
cooperative agreement with the National Science Foundation. 
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In Fig. 2 we show the resulting spectrum. In the analysis of the next section we will also 
make use of some data from the calibration star HD 27697 kindly provided by Peletier et 
al. (1997), obtained in the same spectrograph but with higher spectral resolution. 

3. Results 

Prior to computing the FWHM of the Mg lb lines, the spectrum of the lens galaxy 
had to be continuum normalized and subtracted. First the continuum was drawn using as 
reference the calibration-star spectra and the results of FSMD. We fitted the continuum 
points to splines with CDRAW, a command of DIPSO (Howarth & Murray 1988). Next, 
the spectrum was divided by the fitted continuum. Finally, the continuum level was set to 
zero by subtracting 1. The resulting spectrum (blueshifted dividing A by 1 + z = 1.3569) in 
the wavelength range of the Mg lb triplet is shown in Figure 3. 

To obtain a first estimate of the FWHM we attempted a three- Gaussian fit (constraining 
the central lambdas of the Gaussian functions by using the atomic wavelengths) to the 
Mg lb lines corresponding to Gl. We have supposed the same width for the three lines and 
fitted another two Gaussian functions to the spectral absorption features close to the triplet 
in the red and blue wings, respectively. To restrict the fitting further, we have constrained 
the relative intensities of the lines using as reference values the quotients 0.85:1:0.81, found 
fitting the Mg lb triplet of the calibration star AGK 2+14783. After correcting for the 
instrumental width we obtain a value of 07 = 323 ± 19kms _1 . The uncertainty represents 
the rms error of the multi-Gaussian fit to the data but does not include systematic sources 
of error coming from the use of the reference values for the relative intensity of the lines, 
or from the continuum correction. The last is a very important source of error since the 
galaxy absorption lines are strongly diluted by the QSO emission. 
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To obtain more confident values for 07, we used the standard cross-correlation 
technique, in the following steps: (i) continuum correction (as described above) of the 
galaxy and template spectra, (ii) transformation of the wavelength axis taking logarithms, 
and (iii) cross-correlation of the spectra with an appropriate template using the command 
XCORR of DIPSO, which interpolates the data and masks the edges of each spectrum with 
a cosine bell. As templates we used alternatively: AGK 2+14783 observed on the same 
night as the galaxy under the same instrumental conditions, HD 27697 observed with higher 
resolution, and this last spectrum convolved with a Gaussian of 6.58A (FWHM) to match 
the resolution of the galaxy observations. To derive 07 from the galaxy spectrum/stellar 
template cross-correlatation function, we have used the standard method also followed by 
FSMD. 

This consists in fitting a Gaussian to the central points of the cross-correlation function 
peak. The FWHM of the Gaussian fit depends on the template and galaxy broadenings. 
To find the FWHM (07 ) relationship the method is calibrated by cross-correlating the 
template with itself smoothed with Gaussian functions of different widths. In Fig. 4 we 
present the FWHM(<T;) calibrations corresponding to the three template used by us. Using 
this method we have obtained the following 07 values corresponding to each template: 
SlOi^kms" 1 (template: AGK 2+14783); 310±27km S - 1 (HD27697); and 331 ±35 km s" 1 
(HD27697 convolved with a Gaussian of FWHM = 6.58A ). The uncertainties in the 
previous computations correspond to the rms errors in the cross-correlation fitting width. 
The best estimate corresponds to the star observed on the same night (AGK 2+14783), 
which seems reasonable. Since there are no significant differences with the results obtained 
with the other templates, we adopt it as the final value. In Fig. 5 we present a smoothed 
(a = 310 km s -1 ) spectrum of the calibration star AGK 2+14783 superimposed on the 
0957+561 spectrum. To estimate the global uncertainty, we consider two different sources 
of error: (i) the rms error in the cross-correlation fitting width (17kms _1 ), and (ii) the 
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selection of the template star. The dispersion among the three results derived using the 
three different templates (~ lOkms^ 1 ) is a conservative estimate of (ii). Consequently, we 
derive a total uncertainty of ~ 20kms~ 1 . 

The value for the stellar velocity dispersion of the lens galaxy derived by us (310 ± 20 
km s _1 ) is in agreement with the value of 303 ± 50 km s _1 obtained by Rhee (1991), but 
there is a difference of 31 ± 23 km s _1 between the new result and the averaged value of 
279 ± 12 km s _1 presented by FSMD. To exclude the possibility that this difference arises 
from the data analysis, we repeated the computations following the procedure used by 
FSMD and obtained essentially the same results. There is a relatively good agreement with 
respect to the central stellar velocity dispersion by TF (288 ± 9 km s _1 ). The difference of 
22 ± 22 km s _1 is consistent with zero. 



To relate physically the derived values for 07 (FSMD, TF and this paper) with the 
compact central mass of Gl we focused on a simple scenario for the galaxy: a singular 
isothermal sphere, SIS, with velocity dispersion cr h , plus a central massive dark object, 
MDO, with mass M c . The halo of luminous stars (p{) and the halo of dark objects (p rn ) 
can be considered as spherically symmetric ideal gas spheres (the velocity distribution is 
isotropic). In that case, the gas pressure is 



4. 



Discussion 



4.1. A Simple Isotropic Model 



Pi = Pi(?i, 



(1) 



and the equation of hydrostatic equilibrium, 



dpi/dr = -[GM(r)/r 2 ]p, 



(2) 
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with % = l,m and M(r) = An J r p m r 2 dr. For a SIS model of dark matter, p m = a 2 ./2Gr 2 n 
and <r m = (T/j, a velocity dispersion independent of the radius. However, the inclusion of a 
compact nucleus (SIS+MDO), M(r) = M c + 2a 2 h r/G, leads to a 2 m = tr£ + GM c /3r, and now 
<r m is not constant at every point of the galaxy. To obtain o\ it is necessary to know the 
behavior of pi{r). This information can be inferred by considering that the observed surface 
brightness profile, I(s), traces the surface density of luminous matter, 



Hubble Space Telescope (HST) data indicate that I(ip) oc ip~ n , n « 1 (O'.'l < ip < 5"), 
where the cp — s relationship for Gl (Q = 1) is given by </?(") ^ (l/3)h _1 s(kpc), although 
the whole surface brightness profile is not well fitted by a single power law (Bernstein et 
al. 1997). At the HST resolution limit, the light distribution could be singular, and so, a 
singular isothermal halo (p t = Ar~ 2 ) would be in rough agreement with the HST profile. 
Our master equation will be 



We admit that the previous relations (1-4) are valid at separations (s) and radii (r) much 
greater than the Schwarzschild radius of the mass M c . For example, taking M c = 4xl0 9 h~ l 
Mq, the Schwarzschild radius is of about 4xlCT 4 Yr 1 pc. Therefore, R = 1 h^ 1 pc (ss 1/3 
mas in angular separation) can be a reasonable threshold of validity (we avoid a complex 
description of the region at r < R, which should include relativistic effects, the probable 
existence of an accretion disc, the possible distortion of the luminous isothermal halo as 
due to "cannibalism" by the system central massive dark object + accretion disc, etc.). 

In practice, one does not measure the stellar velocity dispersion, 07 (r), but rather 
either the line-of-sight stellar velocity dispersion at separation s (using a small aperture) or 
the velocity dispersion associated with the luminous stars inside a finite aperture, S. In 




(3) 



of = a 2 h + GM c /3r. 



(4) 
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order to compare with observations, the useful expressions are 



j, m _.W«Ei«<fs 



and 

af(s)=a 2 h + 2GM c /3s7i, (6) 

where £/(s) oc s _1 . If the "forbidden" central circle (irR 2 , where R is the threshold of 
validity; see comment after Eq. (4)) is included in S nR 2 , we assume that its contribution 
to <Ji(S) is negligible. Equations (5-6) are applicable only when there is no seeing, and must 
be properly modified accounting for seeing effects in any realistic case; e.g., 

^ (s ' 5 * } = l^ppa" = aft + (^ao^/ao, (7) 

where P(t,s) oc exp[— (t — s) 2 /2s 2 ] represents the point spread function corresponding to 
a seeing of Gaussian dispersion s* and \ q = J St t~ g P(t,s)d 2 t (q=l,2). Equation (7) is the 
realistic version of Eq. (6). It is now relatively easy to estimate \ q (q=l,2) at s < 2s*, and 
so, to infer 

cr?{s, SJ = a 2 h + y/8/9ir 3 (GM c /s*)(l + In {2s,/ R). (8) 

From this equation, at s < s*, one concludes that the "heating" [AT oc a 2 = crf(s, S*) — a" 2 ] 
caused by the compact nucleus is not traced by a 1/s law [as derived from Eq. (6)]. On 
the contrary, at this region will appear a quasi- isothermal behaviour due to the Gaussian 
smoothing of the seeing. As the angular separations ip are proportional to the separations 
s, we can rewrite a c (in km s _1 ) and 

2 

af(<p, a,) =a 2 h + 2.4 x W~ 7 {hM c / <j*){\ + In (6 x 10V*), (9) 

where hM c is in M , and ip and a, [a, = </?( s *)] are in arcsec. In this paper, we concentrate 
on an isotropic model incorporating a compact central mass, M c , which is responsible for a 
gradient in the stellar velocity-dispersion trend. 



4.2. Application to Real Data 



Using the data derived by FSMD, equation (9) allows us to make a first isotropic 
model measurement of hM c (isotropic model measurements of the mass at the center of the 
galaxies are usual, e.g., Kormendy et al. 1998). We have estimated the angular separation 
(p(i) (i=l,5) corresponding to the centers of the five rows (small rectangles of l"x 072) used 
by FSMD, ranging from </?(3) = 0/21 to </?(5) = 0753. As the spectra were taken with a 
0"8 (FWHM) mean seeing (see Table 1 of FSMD), at ip < 075, the law (9) can be applied 
((j* = 0734). Avoiding the finite size of the rectangles (due to seeing, it is expected that 
the inhomogeneities in the rows will be small), one can establish a relationship between 
the observations 07(2) and [</?(i),<7*] via Eq. (9). Our model fits the five constraints with 
two free parameters, i.e., we have 3 degrees of freedom (dof = 3). The best-fit model gives 
X 2 /dof ~ 3.5, and due to this poor x 2 /dof ', we have drawn the A(x 2 /dof) = 1-10 contours 
rather than those for A% 2 = 1,2,... (see Fig. 6). Our best-fit velocity dispersion is a h = 39 
km s _1 , and models with 07^ > 250 km s _1 correspond to A(x 2 /dof) > 2. This is a first 
surprising result. However, the central MDO has a mass in the interval: (0.5-2.1) x 10 10 hr 1 
M Q [A(x 2 /dof) = 2 limits]. We consider the la, 2a,... bounds such that they correspond 
to Ax 2 = x 2 /dof , Ax 2 = A(x 2 / dof),... rather than Ax 2 = 1,4,... (see, for example, Grogin 
& Narayan 1996a; Barkana et al. 1998), and consequently, here, A(x 2 /dof) = 2 limits are 
1.3a limits (80% confidence limits). So, the FSMD data agree with the existence of a very 
large compact central mass in Gl. 

When the aperture S is large, we must use a modified version of Eq. (5) instead of 
Eq. (7). In the presence of seeing, the velocity dispersion associated with the stars inside a 
finite aperture, S, is given by 



As was shown in §3, we have presented a new measurement of the stellar velocity dispersion. 



af(S, S*) 



f s dSJ St a 2 mi(t)P(t,s)d 2 t 
J s dSJ St E,(f)P(t,s)cPt 



a 2 h + (2GM c /3n)[ 



Is X 2dS 

j s x 1 ds l - 



(10) 
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The mean seeing was 1" (FWHM), and the aperture was a rectangle centered on Gl with 
semisides FWHM/2 = 075 and 3FWHM/2 = l'/5. Under these observational conditions, it 
is possible to show that Eq. (10) can be rewritten as Eq. (9), by using an effective angular 
separation of 2a*/^/3. So, if our measurement is labeled with the number "6", 07(6) will 
be compared to <Ji[<p(6) = 2a*/y/3, a*} = <7j((X'48, 0742). Re-evaluating x 2 /dof with the 
six constraints (dof = 4), we obtain the contour plots appearing in Figure 7. The best-fit 
X 2 /dof, velocity dispersion and central mass are now ps 4, 161 km s _1 and 1.37xl0 10 h~ x 
Mq, respectively Now, adding the constraint inferred from our spectroscopic study, a shift 
of the solutions towards lower masses and larger values of ah appears. 

Finally, we measure hM c using all the dynamical information at ip < 175: 5 constraints 
from FSMD + 1 constraint from this paper + 4 constraints from TF. TF have obtained 
the central stellar velocity dispersion associated with a spectrum which was averaged 
over 11 rows (an aperture of l"x 273) along the slit. Taking into account a point spread 
function characterized by an FWHM m 078 (see TF), the aperture can be approximated 
as a rectangle with semisides FWHM/2 and 3FWHM/2. Then, the averaged measurement 
07(7) must be compared to ai[(p(7) = 2<t*/\/3, a*} = 07 (0739, 0734). However, from the 
results by TF, stellar velocity dispersions of 289.5 ± 7 km s" 1 at <p(8) = O'.'O, 297 ± 10 km 
s- 1 at <p(9) = 076, and 290.5 ± 7 km s" 1 at <p(10) = 077 can be derived. We note that 
Eq. (9) works relatively well even at ip ~ 2c*, and so this analytical law is also applied 
to the last two dynamical data. The new contour plots are shown in Fig. 8, and the new 
best-fit x 2 /dof (dof = 8), velocity dispersion and central mass are 3, 283 km s~ x and 10 9 
hr l M Q , respectively. With the whole data set, M c < 6xl0 9 h' 1 M G at the I.6&7 (90%) 
confidence level, i.e., a result in clear disagreement with the initial one (Fig. 6) based on 
the data of FSMD data. Therefore, we conclude that the existence of a "normal" MDO 
with a mass of about 10 9 h' 1 M Q rather than a very large compact central mass is favored 
by the combined (FSMD + TF + this paper) data set. 
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The predicted values of the stellar velocity dispersions for the final best-fit (solid line) 
and the best-fit to FSMD data (dashed line) are compared to the observational values (la 
limits) in Fig. 9. Filled circles, open circles and an open square represent FSMD data, TF 
data and the new measurement, respectively. The trend observed by TF (open circles) and 
the dynamical value derived in Sect. 3 (the open square) clearly disagree with the solution 
corresponding to the FSMD data (dashed line). From another point of view, three data 
by FSMD are in disagreement with the final fit (solid line) and the data by TF and us. 
Thus, there is evidence in favour of either some kind of error in several measurements or an 
underestimate of several error bars. New studies must throw light on this problem and lead 
to a better reduced chi-squared. 



In the two previous subsections, we assumed an isothermal profile for the luminous 
tracers. However, Fig. 10 gives the observational surface brightness (filled triangles; 
Bernstein et al. 1997) and the results of the profile fitting (solid and dashed lines). If 
I(ip) oc Lp~ n , the data at if < 0'/2 show a slope n = 0.65±0.04 (x 2 /dof m 1), while the 
profile is well fit by a power law with n = 1.26 ±0.01 (x 2 /dof 1) in the interval [0'.'2, 5"]. 
The whole HST trend can be also fitted by a single power law: n = 1.12±0.01. This global 
fit (solid line) is not so well as the partial ones (dashed lines) and suggests the existence of 
a slight global deviation from the isothermal behaviour. 

The switch from n = 1 (isothermal) to n — 1.12 may have some effect on the 
aperture-seeing integrals that appear in the definition of af(S,S*) [see Eq. (10)]. For 1 
< n < 2, one has £;(£) oc t~ n and an interpolated analytical dynamical law 



4.3. 



Some Refinements 



2 




n+1 



(GM c /t). 



(11) 



n + 1 



(n + 2) 2 
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From Eq. (11) and Ej(£) oc t~ n (n = 1.12), and taking the aperture and seeing conditions 
associated with our measurement, we can recalculate the ratio between the two aperture- 
seeing integrals. It is clear that the factors l/(n+l) and (n + l)/(n + 2) 2 in of (£) are almost 
insensitive to the change from n — 1 to n — 1.12 (they vary in 2-6%). On the contrary, an 
accurate numerical estimate indicates a significant change of about 50% in the term related 
to the central mass M c . This last result is due to the use of a brightness profile (Ej oc t~ 1A2 ) 
more realistic than the isothermal profile. 

The high sensitivity of <Ji(S, S*) (in the term proportional to M c ) to the light distribution 
E/, forces us to recompute the central mass. We now use the global value n = 1.12 in 07, the 
double behaviour n = 0.65 at (p < 0'.'2 and n = 1.26 at (p > 0'/2 in Ej, and accurate numerical 
integrals for the ten observational apertures. The new measurements of the central mass 
are in excellent agreement with the "rough measurements" quoted in subsection 4.2. Thus, 
the FSMD data lead to a h < 224 km s" 1 and M c = (3.5l£jj) x lO^/i^M© (lo limits), and 
considering all data, we infer (lo limits) a h = 294^ km s" 1 and M c < 6.2xl0 9 h~ l M & . 
Moreover, the new values of x 2 /dof (best-fits) are very similar to the old ones. 



5. Conclusions 

We have presented new spectroscopy of the lens galaxy (Gl) of QSO 0957+56lA,B and 
have determined the stellar velocity dispersion integrated in a rectangular square of 1" x 3" 
centered on the galaxy. The obtained value, 0/ = 310 ± 20kms~ 1 , is somewhat greater than 
the determination by FSMD (279 ± 12kms _1 ), but in agreement with the measurement of 
288 ± 9kms _1 that TF have published during the preparation of this paper. 

Motivated by the controversy about the amount of compact dark matter at the center 
of Gl (a mass of 4xl0 9 h~ l M & was claimed by FSMD for the central object, whereas TF 
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discounted the presence of a compact nucleus), we have studied the mass of the compact 
nucleus in Gl by means of a simple isotropic model (SIS, singular isothermal sphere, plus 
an MDO) and the available measurements at angular separations < l'/5 (FSMD, TF, and 
this paper). Using exclusively the data of FSMD, one obtains a surprising best-fit mass 
measurement of about 2xl0 10 hr 1 M & for the object in the center of the galaxy. So we find 
that the previous rough estimate by FSMD could differ by almost an order of magnitude 
with respect to the value derived from the SIS+MDO model. At the 80% confidence level 
the compact central mass should be in the range (0.5-2.1) x 10 10 hr 1 M Q . This result would 
be a challenge to current knowledge concerning compact nuclei in galaxies (see, for example, 
fig. 17 of Kormendy et al. 1998). When we consider our measurement in the fit, the 
solution shifts towards lower masses with a best fit of 1.4xl0 10 hr 1 M Q . More significantly, 
adding this new constraint the 80% confidence level includes "normal" MDOs (10 7 -10 9 
M & ) and even the absence of a compact central mass. To complete the analysis, we also 
incorporated the dynamical data by TF. The final fit (using ten different measurements) 
suggests the existence of an MDO with a mass of about (best-fit) 10 9 h~ l M G . At the 
90% confidence level the compact central mass inferred from the ten measurements must 
be < 6xl0 9 h' 1 M & . Finally, we remark that the measurements of the central mass are 
based on the assumption of an isothermal surface density of luminous matter. However, a 
refined model including a more realistic £; (faithfully traced by the observational brightness 
profile) led to very similar results (see subsection 4.3). 

We are especially grateful to Emilio Falco for useful discussions and communications 
and John Tonry for helpful comments about his dynamical results. We also would like to 
thank Gary Bernstein for providing the Hubble Space Telescope surface photometry of the 
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the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de 



-16- 

Canarias. We thank all the staff at the Observatory for their kind support. 

This work was supported in part by DGESIC (Spain) grant PB97-0220-C02. 



REFERENCES 



Barkana, R. et al. 1998, |astro-ph/980809l: 



Bernstein, G., Fischer, P., Tyson, J. A., & Rhee, G. 1997, ApJ, 483, L79 

Bernstein, G., Tyson, J. A., & Kochanek, C. S. 1993, A J, 105, 816 

Borgeest, V., & Refsdal, S. 1984, A&A, 141, 318 

Falco, E. E., Gorenstein, M. V., & Shapiro, I. I. 1991, ApJ, 372, 364 

Falco, E. E., Shapiro, I. I., Moustakas, L. A., & Davis, M. 1997, ApJ, 484, 70 (FSMD) 

Fischer, P., Bernstein, G., Rhee, G., & Tyson, J. A. 1997, AJ, 113, 521 

Greenfield, P. D., Roberts, D. H., & Burke, B. F. 1985, ApJ, 293, 370 

Grogin, N. A., & Narayan, R. 1996a, ApJ, 464, 92 

Grogin, N. A., & Narayan, R. 1996b, ApJ, 473, 570 

Howarth, I. D., & Murray, J. 1988, Starlink User Note 50 

Im, M., Griffiths, R. E., & Ratnatunga, K. U. 1997, ApJ, 475, 457 

Kochanek, C. S. 1991, ApJ, 382, 58 

Kormendy, J., Bender, R., Evans, A. S., & Richstone, D. 1998, AJ, 115, 1823 
Kundic, T. et al. 1997, ApJ, 482, 75 

Oscoz, A., Serra-Ricart, M., Goicoechea, L. J., Buitrago, J., & Mediavilla, E. 1996, ApJ, 
470, L19 

Oscoz, A., Mediavilla, E., Goicoechea, L. J., Serra-Ricart, M., & Buitrago, J. 1997, ApJ, 
479, L89 

Peletier, R. F., Prada, F., Arribas, S., Mediavilla, E., del Burgo, C, Gutierrez, C, 
Garcia-Lorenzo, B., & Vazdekis, A. 1997, Spectrum, 15, 16 



-18- 

Refsdal, S. 1964a, MNRAS, 128, 295 

Refsdal, S. 1964b, MNRAS, 128, 307 

Rhee, G. 1991, Nature, 350, 211 

Schild, R. E. & Thompson, B. 1995, AJ, 109, 1970 

Stockton, A. 1980, ApJ, 242, L141 

Tonry, J. L., & Franx, M. 1998, |astro-ph/9809064] (TF) 



Walsh, D., Carswell, R. F., & Weymann, R. J. 1979, Nature, 279, 381 

Young, P., Gunn, J. E., Kristian, J., Oke, J. B., & Westphal, J. A. 1980, ApJ, 241, 507 



This manuscript was prepared with the AAS IATjrjX macros v4.0. 



Fig. 1. — Slit superimposed on an image of 0957+561 obtained on the same night with the 
auxiliary Cassegrain camera. We have partially removed the A and B QSO images by using 
DAOPHOT. 

Fig. 2. — Spectrum obtained co-adding all the exposures. 

Fig. 3. — Continuum-corrected spectrum in the region of the Mg lb triplet. The spectrum 
has been blueshifted to z = (see text). 

Fig. 4. — Calibration of the FWHM(<T/) relationship for three templates: the continuous, 
dotted, and dashed lines correspond to the AGK 2+14783, HD27697, and HD27697 smoothed 
star spectra, respectively. 

Fig. 5. — Comparison of Gl spectrum with the spectrum corresponding to the calibration 
star AGK 2+14783 smoothed to a = 310km s" 1 

Fig. 6. — Xmm/dof (cross) and A(x 2 /dof) = 1-10 contours (solid lines). We use the five 
dynamical data of Falco et al. (1997). The mass is expressed in terms of M 8 = 10 8 M & . 

Fig. 7. — Xmin/dof (cross) and A(x 2 /dof) = 1-10 contours (solid lines). Our central stellar 
velocity dispersion (six dynamical constraints) has been included. The mass is expressed in 
terms of M 8 = 10 8 M & . 

Fig. 8. — Xmm/dof (cross) and A(x 2 /dof) = 1-10 contours (solid lines). We use all the 
dynamical data (Falco et al. + Tonry & Franx + this paper). The mass is expressed in 
terms of M 8 = 10 s M m . 
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Fig. 9. — The ten measurements (with la error bars) of the stellar velocity dispersion (FSMD 
= filled circles, TF = open circles and this paper = open square) and the predicted values 
in two scenarios: a h = 283 km s" 1 , hM c = 10 9 M (solid line) and a h = 39 km s^ 1 , hM c = 
1.9 x 10 10 M (dashed line). See main text for details. 



Fig. 10. — Surface brightness of the Gl isophotes is plotted vs. the mean radius of the 
isophotes. Filled triangles are the observational values in approximately the V band. Lines 
are different power laws I {(f) oc f~ n : n = 0.65 (dashed line from O'.'l to 0'.'2), n = 1.26 
(dashed line from &!2 to 5") and n = 1.12 (solid line covering the whole range of radii). We 
note that these indexes have been derived from HST observations by Bernstein et al. (1997), 
and the slopes are different to the slope (n ~ 2) inferred from KPNO observations at larger 
radii (Bernstein, Tyson & Kochanek 1993). 
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